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We investigate, via numerical simulations, the phase ordering kinetics of a two-dimensional 
soft-spin 0(2) Ginzburg-Landau model when a reversible mode coupling is included via the 
conserved conjugate momentum of the spin order parameter (the model E). Coarsening of 
the system, when quenched from a disordered state to zero temperature, is observed to be en- 
hanced by the existence of the mode coupling terms. Lhe growth of the characteristic length 
scale L(t) exhibits an effective super-diffusive growth exponent that can be interpreted as 
a positive logarithmic-like correction to a diffusive growth, i.e., L(t) ~ (tint) 1 / 2 . In order 
to understand this behavior, we introduced a simple phenomenological model of coarsening 
based on the annihilation dynamics of a vortex-antivortex pair, incorporating the effect of 
vortex inertia and logarithmically divergent mobility of the vortex. With a suitable choice 
of the parameters, numerical solutions of the simple model can fit the full simulation results 
very adequately. Lhe effective growth exponent in the early time stage is larger due to the 
effect of the vortex inertia, which crosses over into late time stage characterized by posi- 
tive logarithmic correction to a diffusive growth. We also investigated the non-equilibrium 
autocorrelation function from which the so called A exponent can be extracted. We get 
A ~ 1.99(2) which is distinctively larger than the value of A ~ 1.17 for the purely dissipative 
model-A dynamics of non-conserved 0(2) models. 

PACS numbers: 64.60.Ht,64.60.Cn,75.10.Hk,75.40.Gb 



I. INTRODUCTION 



The dynamics of statistical systems quenched from high temperature disordered states to low 
temperature ordered states has been a subject of interest for several decades jl-7]. In typical 
situations, the average length scale L(t) of the ordered domain grows in time as a power law 
L(t) ~ t l l z , where the growth exponent 1/z depends on the dimension of the space and that of the 
relevant order parameters, in addition to the conserved or non-conserved nature of the latter in 

Ft 

the relaxation dynamics [5j. Usually characteristic topological defects such as vortices or domain 
walls are generated in the initial disordered state, and the annihilation of these defects provides the 
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main mechanism of coarsening and phase ordering in the system. The observed self-similarity of 
these coarsening systems at different time instants, is usually represented by the so-called dynamic 
scaling hypothesis of the equal-time spatial correlation function of the order parameter. Dynamic 
scaling hypothesis has been an important ingredient in deriving the properties of phase ordering 
in the late time stage. For example, in combination with the so called energy scaling method 
8], domain growth laws could be extracted for almost all model systems with purely dissipative 
dynamics. 

Previous works on phase ordering focused mostly on the cases where the effect of dissipation 
dominates the coarsening and the motion of the topological defects. On the other hand, in reality, 
there exist various interesting systems exhibiting dynamic processes that cannot be described solely 
by dissipative dynamics. One example can be found in the case of magnetic spin systems, where the 
spins are influenced by neighboring spins via precession interaction terms that are energy conserving 
[9]. These elements in the dynamics are called reversible (non-dissipative) mode coupling. In fact, 
among the various model systems (Model A, B, C, E, F, G, H, J) classified by Hohenberg and 
Halperin[10] that are known to describe the dynamic critical phenomena, only models A, B, and C 
are based on dissipative dynamics alone, with the remaining models E, F, G, H, and J retaining 
reversible mode couplings. Hence, it is quite natural to extend studies on the phase ordering 
dynamics to these model systems. 

Indeed some works have been carried out along this direction. A most familiar example is the 
phase separation dynamics of a binary fluid, in which reversible coupling between the hydrodynamic 
flow and the relative concentration, advection of the order parameter field by the fluid flow, plays 
an important role [lllj Qj- Effects of hydrodynamic flow on the phase ordering kinetics and defect 



dynamics of nematic liquid crystals have been studied 13l4l6fl. Influence of precession on the phase 
ordering of an isotropic Heisenberg magnet in three dimensions (the model J [sj]) has also been 
studied \v\ . 

A particularly interesting example is the phase ordering kinetics of the Bose gas, i.e., time evo- 
lution of the Bose-Einstein condensation, which was studied via the Gross-Pitaevski (GP) equation 
in two and three dimensions, demonstrating the importance of the reversible Josephson precession 
term in the dynamics [is! [lij]. Alternative approach to this problem is to use the appropriate 



stochastic model known as the model F [20j in which the complex order parameter field ip is (both 
statically and dynamically) coupled to the conserved real field m. 

In this work, we investigate the phase ordering process of systems governed by a simpler model, 
namely, the model E [^(J in which there is no static coupling between ip and m. Specifically, we focus 
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our investigation on the effect of the reversible spin precession term on the phase ordering process in 
the soft spin 0(2) models in two dimensions. We especially try to compare the characteristics of the 
phase ordering dynamics of these model systems with those in the case where dissipative dynamics 
alone is considered. The linearized hard spin version of this model was employed by Nelson and 
Fisher 2l( to describe the behavior of the spin wave in the two-dimensional anisotropic ferromagnet 
and the hydrodynamics of the third sound propagation in thin film of He A . In []jJ], the same 
set of linearized equations was used to demonstrate the crucial role played by the non-dissipative 
precession term in the phase ordering process of the defect free case in the two dimensional XY 
model. 

In terms of phase ordering dynamics, 0(2) models in two dimensions are of particular interest 
because the aforementioned energy scaling method do not provide a definitive result on the domain 
growth law. In equilibrium, the two dimensional ferromagnetic XY model exhibits a Berezinskii- 
Kosterlitz-Thouless (BKT) transition at Tbkt due to the unbinding of vortex-antivortex pairs 



22 . 



Below Tbkt, the system has a quasi-ordered phase which is characterized by an algebraic decay of 
the order parameter correlation function for long distances. The critical exponent governing this 
power law decay decreases continuously down to zero temperature. That is, the system is critical at 
equilibrium for all non- vanishing temperatures below Tbkt- Therefore, the coarsening dynamics 
in this model at finite temperatures is expected to exhibit critical dynamic scaling instead of simple 
dynamic scaling. 

Large number of works using the purely dissipative dynamics have been carried out on the 



coarsening dynamics of the two dimensional XY/0(2) models [23M46I] . After all these efforts, it 
is now agreed that, in the phase ordering dynamics of ordinary 0(2) models in two dimensions 
without mode coupling terms and with non-conserved order parameter, the growing length scale 
exhibits a logarithmic correction to the diffusive growth as L(t) ~ (t/ hit) 1 / 2 . Here, the logarithmic 
correction can be attributed to a logarithmic divergence (with the system size) of the effective 
friction constant (or, inverse mobility) of a moving vortex in the dissipative dynamics with non- 



conserved order parameter 
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In a related simulation work [42] on the coarsening dynamics of superconducting Josephson 
junction arrays based on the dynamics of resistively-shunted junction model [52], it was revealed 
that there is no logarithmic correction in the growth law, resulting in purely diffusive growth law. 
This absence of logarithmic correction was understood in terms of finiteness of the effective friction 
constant (or inverse mobility) of a moving vortex in the limit of large system size, which is due to 



the particular type of dissipative coupling in Josephson junction arrays 



42]. 



4 



In addition, the coarsening dynamics of the two dimensional XY model with a purely Hamil- 
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45, 



53H55I] exhibited the growing length scale with an apparent late-time 



tonian dynamics 

power-law growth as L(t) ~ t 1 ^ with the exponent \jz larger than the diffusive exponent 1/2. 
This kind of (dynamics-dependent) non-universal growth law motivated us to investigate further 
the phase ordering dynamics in related model systems. 

Actually, in the model E, we find that the growing length scale L(t) exhibits an apparent super- 
diffusive growth as L(t) ~ t l l z with the "effective" dynamic exponent 1/za little larger than the 
diffusive exponent 1/2. This can be contrasted to the case when reversible mode couplings are 
absent where the apparent growth exponents with values a little smaller than 1/2 are obtained due 
to the (negative) logarithmic corrections mentioned above. The apparent super-diffusive exponent 
in the case of reversible mode coupling may suggest that the asymptotic behavior can be represented 
as a diffusive growth with a positive logarithmic corrections of the form L(t) ~ (tint) 1 / 2 , as 
discussed in detail below. 

One can extract another growing length scale from the excess energy relaxation. It is expected 
that the so-called energy scaling method [8[] can be applied to the present case as well since 
the reversible contribution does not dissipate energy. This method gives the relation AE(t) ~ 
Lg 2 (t) In L/E{t) (between the excess energy AE(t) versus a growing length L%(t)). We find that 
the length scale Le(t) derived from the above relation agrees with the growing length scale of the 
domains L(t) in almost the entire time. 

In order to understand the characteristics of the length scale growth in the present system, 
we noted that the coarsening and the resultant growth of the length scale are determined by the 
annihilation of vort ex- ant i vortex pairs. Therefore, it is plausible to assume that the time t that 
it takes for the system to grow up to a length scale L(t) will correspond to the time that it takes 
for vortex-antivortex pairs of separation L(t) to annihilate. Based on this picture, we devised a 
simple phenomenological model of dynamics for the annihilation of a vortex-antivortex pair where 
we incorporate the effect of vortex inertia and logarithmically divergent mobility of the vortex 



56H58I]. We found that this model can describe closely the simulation results. We argue that the 
microscopic mechanism of these coarsening features including the growth law, is closely related to 
the fact that there exists one conserved quantity in this model. That is, the m field component (see 
the next section), which is conjugate to the planar spin, is conserved. This corresponds to the third 
(z-) component of the spin in Heisenberg spin systems with easy plane anisotropy [21]. Due to the 
conservation of this component, there appear propagating spin wave modes at low temperatures 
below Tkt where the vortices and antivortices are neglected. In the case of coarsening dynamics 
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where there exist a lot of vortices and antivortices in the initial stage, we do not expect that the 
spin waves are fully propagating. We expect rather that, due to the interaction of the vortices 
with the background vortices and antivortices, they propagate only within a short range and then 
scatter off the vortices (and antivortices) in a complicated manner such that the mobility of an 
isolated vortex is enhanced with a logarithmic depence on the size of the system. 

Spin autocorrelation function A(t) also was calculated which is expected to be related to the 
growing length scale L(t) through a new exponent A as A(t) ~ L _A (i). We could extract the 
value of A in the long time limit as A ~ 1.99(2). This result may be interpreted as A = d = 2 
for the present case. . We note that this value is quite distinct from the corresponding value of 
A ~ 1.17 for the case of no reversible mode coupling which was calculated both theoretically as 



well as numerically 
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20, 



3J, 
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60] 



II. THE MODEL E 



We consider the ordering kinetics of the soft spin model known as the model E whose Hamil- 
tonian is given by 



H\6,a] = I d 2 r 



I^f + ^-lf + fm 2 ] (1) 

where (ft is the two-component vector order parameter <p = (0x, <p2) and m is the third component 
conjugate to spin <ft (which corresponds to the z-component of the three-dimensional magnetization 
in Heisenberg spin systems with easy plane anisotropy). When K is set to zero, the above Hamilto- 
nian becomes equivalent to that of 0(2) Ginzburg-Landau model. Therefore, the order parameter 
has 0(2) rotational symmetry on (</>i, </>2)-space and the equilibrium average of m is zero. The 
model- E Hamiltonian, ([1]), does not contain static couplings between the order parameter <p and 
the m field. (On the other hand, the Model F does contain such a static coupling.) However, as 
can be seen below, dynamic couplings between these two fields arise from reversible mode coupling 
terms. The two fields (f and m satisfy the following Poisson bracket relations [l(| 

{<j>i(r),m(r')} = -gM r )$(r - r'), 
{<f>2(r),m(r')} = g<t>i{r)5{r - r'), 

{Mr), Mr')} = (2) 

where g denotes the strength of the mode coupling that is an analog of the gyromagnetic ratio of 
the spins. These Poisson bracket relations generate the reversible mode couplings in the equations 



of motion, which cause spin precession. They play a key role in the critical dynamics [10|]. Note 
that the order parameter eft is not a conserved quantity. Instead, conserved quantity in this system 
is the m component which causes the precession of the order parameter eft in the x — y plane of the 



spin space 



iofi. 



The dynamics of the system described by the fields ip[s (= cfti, <l>%,m) can be written as 

+ Vi(r,t) (3) 



dt ~ ^ 
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where T^j denotes generalized kinetic coefficients (r and D below) and £j denotes the thermal 
noises ( 771, 772 and C below). {A, B} is the Poisson bracket of two generic dynamic variables A and 
B as defined above ([2]). With ([2]) and the above relation, the equations of motion becomes 

dfa SH 6H 

dfo 5H 5H 

W = 9 ^ " r ^ + r?2(r ' t} ' 

where V and D are the kinetic coefficients, and the Gaussian thermal noises ff and £ at temperature 
T satisfy 

(^(r,%(rV)) = 2rr<%<5(r-P)<5(i-i / ), 

(C(r,t)((P,t')) = -2DTV 2 S(f-r')5(t-t'), 

( Vl (r,t)C(r' UO) = (5) 
The equation of motion ([3]) can be explicitly written as 

^ = -gKfom + r[vVi + (i-^i- 0|)0l] +m(r,t), 



ff if0im + r[V 2 2 + (1 - 0? _ ^|)0 2 ] + ^(f, t ), 
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(sTYl — * 

— = gV • (<^iV0 2 - 2 V0i) + D#V 2 m + ((r,t). (6) 

We can see again that, in the limit of g = 0, the order parameter and m are decoupled, resulting in 
the equation of motion for the order parameter equivalent to that of the time-dependent Ginzburg- 
Landau equation for 0(2) model. The equation for m in (|6|) can be written in the form of a 
continuity equation 

dm -> 

~dt = " v ' Jm ' 

J m = <?(</> 2 V0i - 0i V0 2 ) - DKVm + i(r, t) (7) 
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where J m is the corresponding current density. We see from ([7]) that m is a conserved quantity. In 
([7]), the new thermal noise £ satisfies 

c(r,t) = -v-e, 

(&(f, t)^(p, 0) = 2DT6ijS(f-?)6(t - t') (8) 

Note that the equations of motion ([6]) are invariant under the transformation m — > —m, <p\ —> (fe, 
and 4>2 <f>i- 

An important physical element brought on by the reversible mode coupling term is the existence 
of the propagating spin wave especially when we can neglect the existence of vortices and antivor- 



tices. In reference 18|], in the quench dynamics of Gross-Pitaevskii equation, a linear growth of 
z = 1 was found, which was conjectured to be due to the existence of propagating spin waves. Here 
in the case of model E, however, we find only a minor enhancement with logarithmic correction to 
a diffusive growth. This implies that the spin wave here is not fully propagating but interacting in 
some complicated manner with the vortices or antivortices, resulting in suppression of the propa- 
gation due to scattering with the vortices and antivortices. See also the discussions in section IV 
on related questions. 

III. SIMULATION RESULTS AND DISCUSSIONS 

In order to investigate the phase ordering kinetics of the model E, we integrated the spatially 
discretized form of the equations of motion ([5|) at zero temperature T = 0, i.e., in the absence 
of thermal noise, starting from a random initial configuration (corresponding to an equilibrium at 
T = oo). As for the values of the parameters, we put T = D = K = \. In this work, we fix the 
intensity of the precession g = 1. Simulations are performed on systems with spatial discretization 
of square lattice type of dimensions up to 3000 x 3000 with periodic boundary condition. For the 
integration time interval of Euler method, we chose St = 0.05. The simulation results were the 
same for smaller time intervals. 

Two main quantities of interest in phase ordering kinetics of this system are the equal-time 
spatial correlation of the order parameter defined by 

C(r,t) = ^(^2^(t)-^+r(.t)J. (9) 
and the spin autocorrelation function A(t) defined as 
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The equal-time spatial spin correlation function satisfies the simple dynamic scaling (Fig. 1) 

C(r,t) = G(r/L(t)) (11) 

with the correlation length Lit) defined from C(r, i)] r=i (^ = Co- Here Co was conveniently taken 
as Co = 0.4. The correlation length L{t) is found to grow in time as an apparent power law 

L(t)~i* «£~ 0.54(1) (12) 

It is interesting to find that the apparent domain growth exponent in the late time regime 
((f) ~ 0.54) is somewhat larger than the diffusive value of 1/2 (Fig. 2a). This can be contrasted 
to the conventional purely diffusive case (g = 0) where the effective growth exponents obtained 
numerically are invariably smaller than 1/2. Now, in order to understand this behavior, we con- 
struct a simple phenomenological model of coarsening of the system as follows. To begin with, 
we note that the coarsening of the system is dominated by the annihilation of vortex-antivortex 
pairs. It is easy to see that, in order for the system to grow up to a length scale L, vortex pairs 
of sizes on the order of L must be already annihilated. Even though this annihilation process of 
vortex pairs is very complicated many-body process involving many vortex pairs, we assume that 
we can simplify the whole coarsening process (corresponding to the growth of the length scale up 
to L) by the annihilation of single vortex-antivortex pair of size L with suitably defined interaction 
potential and other dynamical parameters. 

From previous works on the dynamics of a vortex in the related model systems of anisotropic 



Heisenberg spin systems (which is based on a collective variable approach) 5614581] , we might assume 
that the vortex acts like a small particle with finite 'mass' moving under the influence of an external 
force with some (length-scale dependent) mobility. Then the distance R between a vortex and an 
anti-vortex would be described by the following equation of motion 

Here we denote the length-dependent effective mass of a vortex as M(R) and the length-dependent 
mobility of a vortex as n{R)- We also assume naturally that the vortex-antivortex pair is interacting 
via a Coulombic force F(R) = —k/R in two dimensions with a proportionality constant k. We will 
set k = 1. 

One important remaining question is how we set the functional form of the M(R) and fJ-(R)- 
From the apparent super diffusive behavior of the growing length scale obtained from the simple 
power law fit with an exponent a little larger than 1/2, we guessed that the growing length scale 
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might be represented, at least in the late time regime, as L(t) ~ (tint) 1 / 2 . If we neglect the mass 
term, this kind of positive logarithmic correction can be obtained by assuming a logarithmically 
divergent vortex mobility fi(R). Similar phenomenon of logarithmically divergent mobility was also 
found in the quasi-two-dimensional diffusion of colloids [61] or diffusion of protein molecules on 



the membranes 62( under hydrodynamic effect, and also in the vortex diffusion in the anisotropic 
Heisenberg system in two dimensions with spin precession As for the form of M(R), we 

turned to the anisotropic Heisenberg model where the effective mass of a vortex with logarithmic 
dependence on the length scale was derived [5614581] . 



Motivated by these phenomena in analogous systems, we assumed that the effective mass and 
the mobility are logarithmically dependent on R as 

M{R) = m + m l hx(R/r ), (i(R) = Ho + m M-R/ro) (14) 

where too, mi, [1q, fM\ are constants and vq denotes a shortest cutoff length scale in the system 
(corresponding to the vortex core size). Note that tjiq corresponds to the effective mass at the 
shortest cutoff length scale and similarly for for the vortex mobility. 

The way we obtain the growth law from this simple model of vortex-antivortex pair dynamics 
is as follows. We first begin with a finite value of the size of vortex-antivortex pair R = R*. For 
convenience, we also set the initial velocity of the vortex to be zero. Then we numerically solve the 
model equation of vortex-antivortex pair dynamics, (13) and (14), to obtain the time r when the 
size of the vortex-antivortex pair R becomes equal to ro the lower cutoff length scale. By plotting 
the resulting relation R* and r we get the growth law of the coarsening dynamics. 

It would be ideal if we could derive the parameters mo, mi, ^o> an d Hi from the full dynamic 
equations, (6) in the absence of thermal noise. But we do not know how to proceed from this 
equation to calculate these parameters. We simply tried to fit our simulation results (especially 
the dependence of the growing length scale on time) with the numerical solutions of our simple 
dynamic model with our tuned values of the parameters such that the two results agree most 
favorably. 

In order to compare the numerical solutions of the simple dynamic model with the simulation 
results, we plot in Fig. 2b L 2 (t)/t versus t from which we can check the behavior of the (multiplica- 
tive) correction to the diffusive growth. It shows that there exists a smooth crossover from early 
time (t < 300.0) behavior with larger slope to later time behavior with smaller (slower) slope. The 
solid line in Fig. 2b represent the results of numerical solutions from the simple model of vortex 
annihilation where we chose the parameters toq = mi = 1.0, /xq = 1-0, and \x\ = 0.225. We see 
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that the numerical solution exhibits a reasonable agreement with the full simulation results (open 
circles). From our numerical solutions to the simple model we see that the early time behavior 
comes from the effect of the inertial term. If the mobility of the vortex is assumed to be a constant 
independent of the size of the vortex-antivortex pair, then the inertial effect ceases to be effective 
in longer time regime and the correction term would approach a constant plateau. However, the 
real simulations of the coarsening exhibits a steady increase in the correction in the late time 
regime. That is why we incorporated a length-dependent mobility of the vortex with a logarithmic 
dependence on the size of the vortex-antivortex pair. We found that, as the coefficient p\ of the 
logarithmic term increases, the later time correction with logarithmic behavior gets stronger (data 
not shown here). 

In the long time limit with large separation between vortices and antivortices, we can expect 
that the inertial effect would be negligible. Therefore, due to the logarithmic divergence of the 
vortex mobility, we can analytically get the dominant asymptotic growth law as 

L{t) ~ (tint) 1 / 2 . (15) 

The numerical solution of (13) in the long time region, shown in Fig. 2b confirms this growth law. 

In addition, Shown in Fig. 3 are the relaxations of the vortex number density p(t) and the 
excess energy density AE(t) in a log-log plot. From the logarithmic slopes, we get apparent 
power-law relaxation in time p{t) ~ i -1 - 03 ! 1 ) and AE(t) ~ t -0 - 95 ^). In order to understand the 
relationship between these relaxation behaviors and also this apparent growth exponent (]12p . we 
note that there exists the so called energy scaling relation between the excess energy AE(t) and 
the length scale L(t) of the domain growth for the case of 0(2) models in two dimensions, i.e., 
AE(t) ~ L~ 2 (t) ln(L(i)/ao) with ao denoting a short distance cutoff corresponding approximately 
to the size of a vortex core. This was first derived by Bray and Rutenberg 8]. Since the reversible 
mode coupling terms satisfy the energy conservation, this relation should be valid in the presence of 
mode coupling as well. From the excess energy relaxation (obtained numerically) we can extract a 
length scale LE(t) through the relation AE(t) ~ L E 2 (t) ln(L^(t)/ao). We can also extract another 
length scale Ly(t) from the defect density relaxation that corresponds to the average separation 
between neighboring defects in such a way as p(t) ~ L v 2 (t) which is based on the assumption that 
the vortices are uniformly (and randomly) distributed in two dimensional space. These two length 
scales as well as the numerical sulution to the aforementioned simple model of vortex-antivortex 
aniihilation, together with L(t) are displayed in Fig. 4a. 

This figure shows that the length scale L e (t) agrees with the growing length scale of the domains 
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L(t) in almost entire times, which implies that the above Bray-Rutenberg relation is valid in 
the zero-temperature coarsening of the model E. The Bray-Rutenberg relation, together with the 
growth law (fl5l) . leads to the asymptotic relaxation behavior for AE(t) as 

^.t-ifl + l^) (16) 



kit 

Figure 4a also shows that there is considerable difference between the length scale Ly(t) ~ p~ 1 ^ 2 (t) 
and the length scale L(t). Although this discrepancy appears to become smaller in the longer 
times, their agreement in the long time limit does not seem to be guaranteed. This feature seems 
to indicate that the assumption of vortices being uniformly distributed is not valid. This is also 
shown in Fig. 4b depicting the vortex number density p(t) versus L(t) which can be fit for a wide 
range of time with p(t) ~ L(t) _1 - 92 W. 

The nonequilibrium spin autocorrelation function A(t) is expected to be related to the growing 
length scale L(t) through a new non-equilibrium exponent A as 

A(t)~L-*(t). (17) 

We could extract the value of A by plotting A(t) versus L(t) as shown in the Fig. 5, where we 
can see that, in the long time limit the value of A approaches A ~ 1.99(2). This value is much 
larger than the value of A ~ 1.17 for the case of non-conserved 0(2) model with no reversible mode 



coupling 



26, 



30, 



33 



59, 



60]. It appears that the higher mobility of the vortices in the present 
model causes a faster decay of the autocorrelation and hence larger value of the A exponent is 
ensued compared to the case of 0(2) model without mode coupling terms. It would be interesting 
to prove the conjecture A = d = 2 analytically. 



IV. SUMMARY AND OUTLOOK 



In this work, we studied the zero-temperature coarsening kinetics of the model E in which the 
order parameter field ip is dynamically coupled with additional conserved field m. We found that 
the phase ordering kinetics of two dimensional 0(2) model is modified considerably due to these 
reversible mode coupling terms. We can summarize the simulation results as follows. The growth of 
typical length scale in the model E is found to exhibit an apparently super-diffusive behavior. We 
introduced a simple phenomenological model for the annihilation dynamics of a vortex-antivortex 
pair incorporating vortex inertia and logarithmically divergent mobility of the vortex. This model 
was shown to describe closely the positive logarithmic corrections of the simulation data, where 
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the inertial effect dominates in the early time stage, while the positive logarithmic effect of the 
vortex mobility dominates in the late time stage. 

We also investigated the autocorrelation function of the order parameter field. The numerical 
result on the value of A exponent for the present case of the model E is approximately A ~ 1.99(2) 
which is quite distinct from the value A ~ 1.17 for the case of purely dissipative dynamics of 0(2) 
model without reversible mode couplings. Even though we do not have an analytic proof, it also 
leads us to conjecture that A = d = 2 for the model E which seems to be closely related to the 
higher mobility of the vortices as compared to the case of no mode coupling. 

We think that the interaction of propagating spin waves (generated by the mode coupling) with 
the vortices (and antivortices) influence the motion of vortices (and antivortices) in such a way 
that the vortex mobility increases logarithmically on the sizes of vortex-antivortex pairs, thereby 
facilitating the annihilation of the vortex-antivortex pairs in the late time stage of coarsening pro- 
cess with a logarithmic correction to a diffusive growth. Further investigation would be necessary 
to understand in more detail the interaction between the vortex defects as well as that between 
the vortices and the propagating spin waves. 

In this work, we only considered the phase ordering kinetics of the model E in two dimensions, 
in which static coupling between the order parameter field and the conserved field m is ignored. 
In view of the recent surge of interest in the phase ordering kinetics in the ultra-cold atomic gases, 
it would be worthwhile to extend our study to the coarsening of the model F which possess a 
static coupling between the order parameter field and the conserved third component. In Ref. \)\ . 
the superfluid ordering kinetics of the Bose gas followed by an instantaneous quench from high- 
temperature normal state was studied using the Gross-Pitaevskii equation for the order parameter 
field in both two and three dimensions. The order-parameter correlation function was shown to 
obey a critical dynamic scaling for d = 2 and a simple scaling for d = 3 with the growing length 
scale exhibiting a linear growth in time, i.e., L(t) ~ t l l z with the dynamic exponent z ~ 1.0 due 
to the presence of propagating spin waves in both dimensions. Therefore it would be interesting 
to carry out a detailed study on the coarsening kinetics of the corresponding stochastic model, the 
model F, in order to see whether the observed linear growth in the Gross-Pitaevskii equation is 
preserved in the model F or not. 

Recently, the Bose-Einstein condensate with nonzero spin, i.e., the spinor condensates have 



been experimentally realized for alkali atoms such as Na 



63 1 and 87 Rb 64(]. Interplay between 



superfluid and magnetic ordering in these spinor condensates lead to unusual phases and the cor- 



responding topological defects 



651 ] . One remarkable example is the existence of the nematic super- 
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fluid phase and the associated half- vortices in the so-called polar phase (observed in sodium atoms), 



which are usually observed in nematic phase of liquid crystals [66j, [67| . Quite recently, the model F 
dynamics was generalized and applied to the magnetic domain growth in the ferromagnetic phase 
of spinor condensates (observed in rubidium atoms), with and without the constraint of the total 



magnetization conservation 



68 



69J]. Extending our study along these fascinating developments in 



the nonequilibrium dynamics of ultracold atoms would be quite rewarding. 
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(a) 
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FIG. 1: (a) The equal time spatial correlation functions for the 0(2) order parameter of the model E at 
various time stages (t = 10, 20,40,80, 160,320,640,1280,2560,5120) with the system size 3000 x 3000 and (b) 
the scaling collapse of the data in (a) with the appropriate scaling length L(t) for t = 320,640,1280,2560, 
and 5120 as explained in the main text. 




FIG. 2: (a) Growing length L(t) vs. t (open circles) manifesting an apparent power law growth L(t) ~ 
with the effective growth exponent of <f> ~ 0.54(1) (dotted line), (b) L(t) 2 /t vs. t from the full simulation 
(open circles) is displayed here together with the result of a numerical solution (solid line) of simple dynamic 
model of vortex-antivortex pair annihilation with a suitably chosen set of parameters that agrees reasonably 
well with the simulation results. Note that the result from the simple model is scaled by a constant factor 
to fit the full simulation results. 
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10° 




FIG. 3: Relaxation of the vortex number density p(t) (solid line) and the excess energy AE(t) (dashed line) 
exhibiting p(t) ~ f" 10 ^ 1 ) and AE(t) ~ i- - 9 ^ 1 ) respectively 




FIG. 4: (a) The growing length L(t) (open circles) is shown together with the two length scales Ls(t) (dotted 
line) , Ly (t) (dashed line) derived from the energy scaling relation and the decay of the vortex number density, 
respectively, in addition to the numerical solution R(t) (solid line) of the vortex annihilation model as shown 
in Fig. 2b. We find that i_e(i) and R(t) agree very closely with L{t). (b) The vortex number density p(t) 
versus L(t) exhibits a power law relation with p(t) <~ ^(t) -1 - 92 ' 1 ) consistent with the small discrepancy 
between L(t) and Ly{t) in (a). 
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FIG. 5: Nonequilibrium autocorrelation function A{t) of the 0(2) order parameter versus the growing 
length scale L(t) (open circles), which exhibits a relation of A(t) ~ L(t) x with A ~ 1.99(2) in the long time 
limit. 



